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Abstract 

Soliton dynamics in a large variety of longitudinally modulated lattices are studied in terms of 
phase space analysis for an effective particle approach and direct numerical simulations. Complex 
soliton dynamics are shown to depend strongly on both their power/width and their initial mo- 
mentum as well as on lattice parameters. A rish set of qualitatively distinct dynamical features of 
soliton propagation that have no counterpart in longitudinally uniform lattices is illustrated. This 
set includes cases of enhanced soliton mobility, dynamical switching, extended trapping in sev- 
eral transverse lattice periods, and quasiperiodic trapping, which are promising for soliton control 
applications. 
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I. INTRODUCTION 



The propagation of spatially localized waves in inhomogeneous nonlinear media is a sub- 
ject of general interest that rises in many branches of physics such as light propagation 
in photonic lattices, matter-wave formation in Bose-Einstein condensates and solid state 
physics. Periodic modulation of the material properties in the form of a photonic lattice 
configuration has been shown to alter solitary wave formation and dynamics so that lattice 
solitons have features and properties with no counterpart in homogeneous media The 
inhomogeneity of the medium can be either written in the form of waveguide arrays and 
optical networks or dynamically formed and leads to additional functionality with respect 
to wave control allowing for a variety of interesting applications. In such lattices stable 
solitons can be formed in specific positions determined by the lattice geometry. For the 
case of photonic structures transversely modulated by monochromatic (single wavenumber) 
variations of the linear or nonlinear refractive index stable solitons are always formed in 

nn 

the positions corresponding to the minima of the respective potential [2|, |3[. In the case of 
polychromatic lattices the effective potential is different for solitons having different power 
and spatial width so that the number and the positions of stable solitons depend on their 
properties Therefore, the increased complexity of the medium modulation results in 
additional functionality of the photonic structure with respect to soliton discrimination. 

An additional degree of freedom in lattice modulation is related to the longitudinal pe- 
riodic modulation of the properties of the medium. Such cases have been studied in the 
motion of charged-particle wave packets along periodic potentials in an a.c. electric field 
^-[6| and in matter-wave dynamics in Bose-Einstein condensates under the influence of a 
time-periodic potential [t- 9] . In the context of light propagation in photonic lattices, pre- 
vious studies include diffraction managed solitons 10| as well as Rabi-like oscillations and 
stimulated mode transitions with linear and nonlinear wave states in properly modulated 
waveguides and lattices [ll|. Longitudinally modulated lattices have been considered for 
soliton steering in dynamically induced lattices 12Hl4j|. Dynamical localization and control 
has also been studied in longitudinally modulated waveguide arrays including periodically 



curved [15( waveguides and waveguides with periodically varying width [16| as well as lattices 
with longitudinally varying refractive index 

In this work we investigate the effect of the a wide variety of different types of longitudinal 
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lattice modulation on soliton dynamics and explore the additional functionality of the cor- 
responding photonic structure in terms of soliton control based on both their power /width 
and momentum. Strong momentum dependence of soliton dynamics results from the fact 
that longitudinally modulated lattices actually carry momentum if seen as wave modulations 
of the medium [l2j. The two types of longitudinal modulations considered are amplitude 
and transverse wavenumber modulation along the propagation distance, being capable of 
describing several realistic configurations as well as providing fundamental understanding 
of the dynamical features related to even more general modulations. Soliton propagation 
in such inhomogeneous media is modeled by a NonLinear Schrodinger (NLS) equation with 
transversely and longitudinally varying linear refractive index. Utilization of this continu- 
ous model allows for the study of soliton dynamical trapping and detrapping dynamics in 
contrast to discrete models where discrete solitons are locked at high powers to their input 
waveguides and are not allowed to travel sideways. 

In contrast to many previous studies where longitudinal and transverse length scales are 
well separated, our study focuses on cases where length scale interplay takes place. There- 
fore, the periods of longitudinal lattice modulation are comparable to periods of soliton 
oscillation in a longitudinally homogeneous lattices. Resonances between these periodicities 
result in drastic modification of soliton dynamics and gives rise to a plethora of novel dy- 
namical features depending on soliton power and momentum, including enhanced mobility, 
dynamical switching and trapping in extended areas of the lattice, periodic and quasiperiodic 
oscillations. An effective particle approach 18[ is utilized in order to obtain a nonintegrable 
Hamiltonian system describing soliton " center of mass" motion. The complex dynamics of 
the system are analyzed in terms of Poincare surfaces of section providing a comprehensive 
illustration of soliton dynamical features. The remarkable agreement of direct numerical 
simulations with the phase space analysis of the effective particle dynamics suggest that 
the latter provides a useful tool for analyzing dynamics of different solitons propagating in 
lattices with different configurations as well as for designing lattices with desirable features. 
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II. MODEL AND EFFECTIVE PARTICLE APPROACH 



The wave propagation in an inhomogeneous medium with Kerr-type nonlinearity is de- 
scribed by the perturbed NonLinear Schrodinger (NLS) equation 

du d 2 u 

2— — h — — - + 2|upu -|- cn(x, z)u = (1) 
oz ox 2 

where z and x are the normalized propagation distance and transverse coordinates respec- 
tively. We consider relatively weakly inhomogeneous media where the potential function 
n(x, z) models longitudinally modulated lattices, e is a small dimensionless parameter indi- 
cating the strength of the potential and n(x, z) is periodic in x and z. 

The lattice potential has a fundamental transverse (x) periodicity (Kq 1 ) and we consider 
two types of longitudinal (z) modulations. The first one will be called AM modulation, as 
we modulate the amplitude of the potential function according to 

tiam(z, z) = [A + a sin(J^ix + 0) sin(f2z)] sin(K x) (2) 

In Figs. l(a)-(c) we demonstrate three examples of such lattice patterns in the (x, z) plane. 
For the second type of modulation, which we call WM, the transverse wavenumber of the 
lattice varies periodically with z as follows 

%jf(x, z) = Aq sinfi^o^ + ot sin(i^ix + 0) sin(fi^)] (3) 

and the corresponding patterns are shown in Figs. l(d)-(f). Aq is the strength of the 
unperturbed transverse lattice which is set to unity without loss of generality, and a can 
be seen as the relative strength of the AM or WM modulation. The wavenumber of the 
longitudinal modulation is Q and the secondary transverse wavenumber K\ corresponds to 
a nonuniform longitudinal modulation depending also on the transverse coordinate, thus 
allowing for the description of a very wide variety of lattice patterns: We are able to model 
modulations where the positions of the potential maxima and minima are either varying 
in-phase or out-of-phase as in Figs. l(a)-(c) or are curved as in Figs. 1 (d)-(f ), with respect 
to z. 

The unperturbed NLS (e = 0) has a fundamental singe-soliton solution of the form 



u(x, z) = r]sech[r](x — xo)}e l(2 
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where xq = v, a = —v 2 /8 + r] 2 /2. 77 is the amplitude or the inverse width of the soliton 
solution, xq its center (called sometimes the center of mass due to the effective-particle 
analogy of the solitons), v the velocity, and a the nonlinear phase shift. 

The longitudinal evolution of the center x under the lattice perturbation is obtained 
by applying the effective-particle method. According to this method we assume that the 
functional form and the properties (width, power) of the soliton are conserved in the case of 
the weakly perturbed NLS. This assumption has to be verified through numerical integration 
of Eq. ([!]) at least to a good approximation. The equation for x is equivalent to an equation 
that describes the motion of a particle under the influence of an effective potential (periodic 
in our case) and is given by 

dV eff , A . 

mxQ = (4) 

ox Q 

where m — J \u\ 2 dx the integral of the non-dimensional soliton power, equivalent to the 
particle mass, and V e ff(x ) = 2 J n(x, z)\u(x; x )\ 2 dx the aforementioned effective potential 
determining soliton dynamics. By substituting the expressions of fundamental soliton and 
lattice potential Eq. (J2J) or Eq. (j3J), we obtain m = 2rj and 

r (AM) _ nA _^K sm(K x ) 



f^ ± sinh(7riV2^) 1 1 



with K± — K\ ± K , and 



e// )=2 ^E4(^mN) sinh (7 riW2,) () 

with K m = mK\ + K , for the two modulation types respectively. 

The soliton dynamics, as described by Eq. (jlj), with the above effective potentials are 
determined by a nonautonomous Hamiltonian system 

H = ^f + V eff (x,z) (7) 

where z is considered as "time" and v = xq is the velocity of soliton center of mass. 

For the case of zero longitudinal modulation (a = 0) the system is integrable and soliton 
dynamics are completely determined by the form of the z- independent effective potential. 
In such case solitons are either trapped, oscillating between two maxima of the effective 
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potential, or traveling along the transversely inhomogeneous medium. For trapped solitons 
the oscillation frequency varies from a maximum frequency 



corresponding to small harmonic oscillations around the minimum of the effective potential 
to a minimum zero frequency (infinite period) corresponding to an heteroclinic orbit con- 
necting the unstable saddle points located at the maxima of the effective potential. The 
heteroclinic orbit is the separatrix between trapped and traveling solitons. Conditions for 
soliton trapping are determined by the initial soliton energy H, which depends on both soli- 
ton position and velocity, and occurs when —Uq/Kq < 2t]H < Uq/Kq. Moreover, solitons 
located at the stable points (corresponding to the minima of the effective potential) can be 
detrapped if their initial velocities exceed a critical value 



and travel at a direction determined by the sign of their initial velocity. 

The presence of an explicit z dependence in the effective potential [a ^ 0) results in the 
nonintegrability of the Hamiltonian system which describes the soliton motion and allows for 
a plethora of qualitatively different soliton evolution scenarios. The corresponding richness 
and complexity of soliton dynamics opens a large range of possibilities for interesting appli- 
cations where the underlying inhomogeneity results in advanced functionality of the medium. 
The nonintegrability results in the destruction of the heteroclinic orbit (separatrix), allow- 
ing for dynamical trapping and detrapping of solitons. Therefore, we can have conditions 
for enhanced soliton mobility: Solitons with small initial energy can travel through the lat- 
tice as well as hop between adjacent potential wells and dynamically be trapped in a much 
wider area, including several potential minima. Moreover, resonances between the frequen- 
cies of the internal unperturbed soliton motion and the ^-modulation frequencies result in 
quasiperiodic trapping and symmetry breaking with respect to the velocity sign. It is worth 
mentioning that all these properties depend strongly on soliton characteristics, namely 77, so 
that different solitons undergo qualitatively and quantitatively distinct dynamical evolution 
in the same inhomogeneous medium. 




nK /2r] 



(8) 





(9) 
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III. RESULTS AND DISCUSSION 



In the following, we consider relatively weakly modulated lattices being of interest in 
most applications and set e = 10 -2 . For example, in typical optical lattices consisting of 
nonlinear material of AlGaAs type (refractive index n ~ 3.32 at Ao — 1.53/im), when the 
trasnverse dimension x is normalized to l-i-3/im, e = 10 -2 corresponds to an actual refractive 
index contrast An = 10 -4 -j- 10~ 5 which is relevant for experimental configurations. For 
such normalized propagation distance z = 500, which is large enough in order to 

observe the dynamical features presented in the following, corresponds to an actual length 
of 13 -T- 123mm. Moreover, as we show in following, weak lattices are characterized by 
complex but not completely chaotic soliton dynamics, allowing for controllable evolution 
features. It is expected that interesting soliton dynamics occur when the spatial scales 
of the system, namely transverse and longitudinal modulation periods and soliton width, 
become comparable, so that length scale competition takes place. As shown from Eq. (jSj), 
the maximum oscillation frequency of a trapped soliton in an unmodulated lattice depends 
strongly on the ratio K /2rj, so that solitons with different rj can have different ujq as depicted 
in Fig. 2, where it is shown that the value of u saturates quickly to an upper bound 
u o{v oo) = y/2eA K as the width (t/" 1 ) of the soliton becomes comparable or smaller 
than the fundamental transverse period (2ir/K ) of the unmodulated lattice. The value of 
u)q for each soliton and K Q is crucial since by introducing the longitudinal modulation, we 
expect the frequencies Q that have significant effect in the soliton behavior are those that 
fulfill a resonance condition with the unperturbed frequencies (a = 0) ranging from to 
ujo- These resonant interactions give rise to qualitative and quantitative different evolution 
of (xq,v). In order to investigate the new features the longitudinal modulation induces, we 
use the phase space representation of (jl]) on a Poincare surface of section stroboscopically 
produced with period Q. The phase space topology visualizes all the qualitative features of 
soliton dynamics and provides a significant amount of comprehensive information which is 
useful for categorizing and analyzing characteristic cases and conceptual design of potential 
applications. 

As shown in Figs. 3 and 4, in both AM and WM lattices the phase space topology changes 
significantly from the corresponding case with no longitudinal modulation (depicted with 
solid curves). The topology depends strongly on both the characteristics of the lattice and 
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the characteristics of the soliton (77), so that not only different lattices but also different 
solitons on the same lattice can have drastically different dynamical features. This is shown 
in different columns and rows of Figs. 3 and 4, for transverse lattice modulations of period 
A = 27t(A' = f) and soliton parameter 77 = 0.5 ~ 2.5 corresponding to a FWHM value 
range 5.3 -j- f in normalized transverse dimensions. Before proceeding to the investigation 
of specific characteristic cases of soliton evolution we discuss the qualitative topological 
characteristics of phase spaces corresponding to different lattices and soliton power (77). As 
expected for a nonintegrable system, the typical Poincare surface of section consists of reg- 
ular curves corresponding to perturbed tori of nonresonant quasiperiodic orbits, resonant 
islands around periodic closed orbits and densely filled chaotic regions related to complex 
motion. In all cases, the separatrix between bounded and unbounded motion has been 
destroyed and replaced by a chaotic region while the regions corresponding to bounded mo- 
tion in longitudinally modulated lattices are downsized. The extent of these areas depends 
strongly on the amplitude of the effective potential which varies exponentially on the ratio 
K m /r) as shown in Eqs. (jSJ) and fl6]). Therefore, solitons having different power (77) have 
completely different trapping conditions in the same lattice, while soliton mobility can be 
drastically enhanced. In addition, changing the soliton power for a specific lattice config- 
uration results in a drastic change on the frequency spectrum of the unperturbed soliton 
oscillations through (jJq{t]). This leads to the possibility of fulfilling resonance conditions with 
the frequency of the longitudinal lattice modulation Q and to the appearance of periodic or- 
bits and surrounding resonant islands. Another qualitative differentiation between the phase 
spaces is the bifurcation of the stable center corresponding to the minimum of the effective 
potential in the unmodulated lattice to a saddle point seen, for example, in Fig. 3(c),(d) 
for 77 = 1.75,2.5 (Q/ujq = 1,0.97 correspondingly). In addition to the power dependence 
of the phase space structure, a common characteristic of all phase spaces for longitudinally 
modulated lattices is the symmetry breaking of the Poincare surface of section with respect 
to v — and KqXq = —n/2. This feature reveals a selectivity of the lattice on initial velocity 
(momentum) and displacement: Solitons with opposite velocities or symmetrically placed 
from the former minimum of the effective potential (which in the unmodulated case would 
remain in the same orbit) they undergo now qualitatively distinct evolution. In such cases 
one of the solitons can be trapped while the other is detrapped. The velocity selectivity 
is a direct consequence of the momentum that is incorporated in the lattice pattern due 
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to the biperiodic lattice potential and even though is met in both AM and WM lattices it 
seems to be more prominent in the later as seen from the Poincare surfaces of section. The 
dependence of motion on its initial displacement is related to the fact that the local minima 
of the transverse lattice profile changes periodically with z. 

Having discussed the topological features of phase spaces corresponding to different com- 
binations of lattice configurations and soliton powers (rj), we now show specific characteristic 
cases of soliton motion having qualitatively distinct properties and being promising for po- 
tential applications. Although similar cases can be met in many different cases, we focus our 
analysis on the cases depicted in Figs. 3(c) and 4(h) with parts of them illustrated in more 
detail in Figs. 5(a) and (b). In the following figures white thick lines and black dashed lines 
depict soliton center motion as obtained from numerical simulation of the original perturbed 
NLS Eq. ([1]) and the effective particle approach respectively, showing a remarkable agree- 
ment. The thick solid lines depict numerical simulations for the case of an unmodulated 
lattice for comparison. 

In Fig. 6(a) we illustrate soliton propagation for an initial center position and velocity 
corresponding to point (i) of Fig. 5(a). It is shown that soliton undergoes a complex 
evolution being dynamically trapped between several transverse lattice periods, in contrast 
to the case of same initial conditions in an unmodulated lattice. Soliton mobility is even 
more pronounced in the case shown in Fig. 6(b) where complete dynamical detrapping 
takes place allowing for a soliton to travel across the lattice. The latter is the typical case 
of soliton evolution for initial conditions located outside the region occupied by regular 
orbits in the corresponding phase spaces. It is worth mentioning that the degree of mobility 
enhancement is not the same for different solitons in the same lattice as it can be seen by 
the comparison of the size of the regular area in Figs. 3(a)-(d). Initial conditions located 
on the regular or island curves of the phase space, as point (iii) in Fig. 5(a), correspond to 
quasiperiodic soliton oscillations as shown in Fig. 6(c). In Fig. 6(d) the soliton evolution 
for an initial condition close to the saddle point (iv) of Fig. 5(a) is depicted exhibiting an 
unstable (hyperbolic) type of periodic orbit. 

A very interesting evolution scenario is illustrated in Fig. 7(a) for initial conditions 
depicted by points (i) and (ii) of Fig. 5(b), corresponding to solitons having the same ini- 
tial positions but opposite initial velocities: The soliton with the positive velocity remains 
trapped and periodically oscillating (since it corresponds to a center of a resonant island - ex- 
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act resonance) while the soliton with the negative velocity undergoes a dynamical switching 
to a neighbor lattice position where undergoes persistent trapping. This type of feature can 
be considered for promising power and velocity dependent soliton switching applications. In 
the same fashion, the effect of symmetry breaking with respect to v = can lead to velocity 
sign dependent soliton trapping or traveling across the lattice, as shown in Fig. 7(b). A 
trapped soliton propagation having the form of a beat is depicted in Fig. 7(c) for an initial 
condition corresponding to a resonant island. 

As a final case we consider soliton dynamical trapping within an extended area including 
many transverse periods of the lattice in a persistent periodic fashion, as shown in Figs. 
8(b) and (c). The initial conditions leading to such evolution correspond to two families of 
interconnected resonant islands (1) and (2) shown in Fig. 8(a). These islands are located 
outside the separatrix of the respective unmodulated lattice showing that outside but close 
to the separatrix the longitudinal modulation can induce interesting persistent dynamical 
trapping for initial conditions for which traveling solitons are expected in the unmodulated 
lattice. In accordance with previous cases, the lattice possesses a selectivity property with 
respect to initial soliton velocity direction as we have not a symmetry with respect to the 
initial velocity axis. The type of evolution depicted in Figs. 8(b) and (c) is interesting for 
power and momentum dependent multi-port soliton switching applications, where different 
input/output ports correspond to different potential minima. 

IV. CONCLUSIONS 

We have studied soliton dynamics in a large variety of longitudinally modulated lattices 
in terms of direct numerical simulations as well as phase space analysis for an effective 
particle approach. The remarkable agreement of the results suggest that the effective particle 
approach and the phase space analysis with the utilization of Poincare surfaces of section 
provides a useful tool for studying complex soliton dynamics in such lattices as well as 
analyzing and/or designing lattices having desirable properties. It is shown that soliton 
dynamics depend strongly on their power and the corresponding spatial width through its 
relation with the transverse lattice period as well as on both magnitude and direction of 
their initial velocity. A large variety of qualitatively distinct dynamical features of soliton 
propagation have been shown that have no counterpart in longitudinally uniform lattices. 
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Therefore, cases of enhanced soliton mobility, dynamical switching and trapping in several 
transverse lattice periods as well as quasiperiodic and periodic trapping have been shown, 
suggesting that the corresponding complexity of the effective particle phase space gives rise 
to a plethora of dynamical features which are promising for applications. 
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FIG. 1: (Color online) Lattice patterns in (x,z) plane for AM (a-c) and WM (d-f) modulation. 
K = 1 and Q = 0.1 in all cases and K\ = 0, <f) = vr/2 (a,d), K 1 = l,4> = (b,e), K x = 1/2, <j) = 
(c,f). Blue (light gray) colored areas indicate negative values and contain the minima of the lattice 
potential while red (colored areas indicate positive values where the corresponding maxima are 
located. 
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FIG. 2: The spatial frequency of trapped soliton small oscillations ojq as a function of inverse width 
rj for different values of the transverse wavenumber Kq in an unmodulated lattice. 
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FIG. 3: (Color online) Poincare surfaces of section of the effective-particle system for AM lattice 
modulations. Each row corresponds to the same lattice pattern (Figs.l(a)-(c)) and each column 
to the same soliton width with values r/ = 0.5, 1, 1.75, 2 from left to right . In all cases Kq = 1. 
From top row to bottom, we have O = 0.1325, f2 = 0.4, Q = 0.15 correspondingly. All surfaces 
of section are superimposed on the corresponding phase space without longitudinal modulation 
(continuous, red curves). The separatrix between the trapped and untrapped soliton motion in the 
unmodulated lattice (blue curve) is also shown. 



15 




FIG. 4: (Color online) Poincare surfaces of sections of the effective-particle system for WM lattice 
modulations. Each row corresponds to the same lattice pattern (Figs.l(d)-(f)) and each column to 
the same soliton width with values rj = 0.5, 1, 1.75, 2 from left to right . From top row to bottom, 
we have = 0.1, Q = 0.4, f2 = 0.175 correspondingly. All surfaces of section are superimposed 
on the corresponding phase space without longitudinal modulation (continuous, red curves). The 
separatrix between the trapped and untrapped soliton motion in the unmodulated lattice (blue 
curve) is also shown. 
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FIG. 5: (Color online) (a) Detail of the Poincare surface of section from Fig. 3(c) (AM modulation). 
The values of the initial parameters (xq,u) at the points depicted in figure are (i)(—7r/2, 0.202), 
(m)(-1. 57, -0.1425), (ra)(-l,-0.21), (iv)(-l. 5308, -0.005). (b) Detail of the Poincare surface of 
section from Fig. 4(h), (WM modulation). The values of the initial parameters (xo,u) at the points 
depicted in figure are (t), (**)(— 2.5, ±0.07), (Hi), (tu)(-0.34, ±0.085), (v) (-1.58, 0.046). 
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FIG. 6: (Color online) Evolution of the soliton center in an AM lattice pattern with initial condi- 
tions taken from the corresponding points (i) — (iv) of Fig. 5(a). (a) point (i): enhanced mobility 
with dynamic trapping, (b) point (ii): enhanced mobility with complete detrapping (c) point {in): 
quasiperiodic oscillations and (d) point (iv): hyperbolic periodic oscillations. Thick white curves 
correspond to results from direct numerical simulations, dashed black curves from the effective 
particle model and solid black ones from direct simulations for the unmodulated lattice. The 
corresponding lattice pattern is shown in the background. 
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FIG. 7: (Color online) Evolution of the soliton center in an WM lattice pattern with initial con- 
ditions taken from the corresponding points (i) — (v) of Fig. 5(b). (a) point (i) (upper, yellow 
curve): Trapping with periodic oscillations for the soliton with the positive initial velocity, point 
(ii) (lower, white curve): Dynamic switching for the soliton with the negative velocity, (b) point 
(in) (upper, yellow curve): Trapping with quasiperiodic oscillations for the soliton with positive 
initial velocity, point (iv) (lower, white curve): Detrapped motion for the soliton with the negative 
velocity, (c) point (v) (white curve): Beat oscillations. Thick white/yellow curves correspond to 
results from direct numerical simulations, dashed black curves from the effective particle model 
and solid black ones from direct simulations for the unmodulated lattice. The corresponding lattice 
pattern is shown in the background. 
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FIG. 8: (Color online) (a) Poincare surface of section for a soliton with 7] = 5 at a WM lattice 
potential of Fig. 1(d), superimposed on the phase space of the unmodulated lattice. For illustration 
purposes, the surface of section is produced with initial values at one of the four periods of the 
phase space. (1) and (2) are two different families (belonging to different tori) of interconnected 
resonant islands, (b) Soliton center trajectories of each fammily with initial conditions (1) (xq, v) = 
(0.05,0.178) and (2) (xo,v) = (— vr/2, — 0.28). Thick white curves correspond to results obtained 
from direct numerical simulations, dashed black curves are obtained form the effective particle 
model. 
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